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Abstract 



The discrepancy in the past years of many more black-hole soft X-ray transients 
(SXTs), of which a dozen have now been identified, had challenged accepted wisdom 
^ \ in black hole evolution. Reconstruction in the literature of high-mass X-ray binaries 

■ has required stars of up to ~ 40 Mq to evolve into low-mass compact objects, setting 

this mass as the limit often used for black hole formation in population syntheses. 
^ \ On the other hand, the sheer number of inferred SXTs requires that many, if not 

' most, stars of ZAMS masses 20 — 35 end up as black holes (Portegies Zwart et 

al. 1997; Ergma & van den Heuvel 1998). 

In this paper we show that this can be understood by challenging the accepted 
wisdom that the result of helium core burning in a massive star is independent of 
whether the core is covered by a hydrogen envelope, or "naked" while it burns. The 
latter case occurs in binaries when the envelope of the more massive star is trans- 
ferred to the companion by Roche Lobe overflow while in either main sequence or red 
giant stage. For solar metallicity, whereas the helium cores which burn while naked 
essentially never go into high-mass black holes, those that burn while clothed do so, 
beginning at ZAMS mass ~ 20 Mq, the precise mass depending on the ^^C(a, 7)^^0 
rate as we outline. In this way the SXTs can be evolved, provided that the H enve- 
lope of the massive star is removed only following the He core burning. 

Whereas this scenario was already outlined in 1998 by Brown, Lee, &: Bethe (1999) 
their work was based on evolutionary calculations of Woosley, Langer, & Weaver 
(1995) which employed wind loss rates which were too high. In this article we collect 
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results for lower, more correct wind loss rates, finding that these change the results 
only little. 

We go into the details of carbon burning in order to reconstruct why the low Fe 
core masses from naked He stars are relatively insensitive to wind loss rate. The 
main reason is that without the helium produced by burning the hydrogen envelope, 
which is convected to the carbon in a clothed star, a central ^^C abundance of ~ 1/3 
remains unburned in a naked star following He core burning. The later convective 
burning through -"^^C+^^C reactions occurs at a temperature T ~ 80 keV. 

Finally, we show that in order to evolve a black hole of mass ^ 10 such 
as observed in Cyg X-1, even employing extremely massive progenitors of ZAMS 
mass ^ 60 Mq for the black hole, the core must be covered by hydrogen during a 
substantial fraction of the core burning. In other words, the progenitor must be a 
WNL star. We evolve Cyg X-1 in an analogous way to which the SXTs are evolved, 
the difference being that the companion in Cyg X-1 is more massive than those in 
the SXTs, so that Cyg X-1 shines continuously. 
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1 Introduction 



Accepted wisdom in the astropliysical community (for recent relevant articles 
see Portegies Zwart, Verbunt, & Ergma 1997; Ergma & van den Heuvel 1998) 
has been that the results of burning the helium core of a star arc independent 
of whether it is covered by a hydrogen envelope or whether it burns as a 
"naked" star, the envelope having been removed; e.g. by Roche Lobe overflow 
in main sequence or red giant stage of the star. This latter situation takes place 
in binary evolution where the evolution referred to in the above two articles 
determines the black hole progenitors to have ZAMS mass > 40 M©. Thus, 
this mass limit for black hole formation is adopted for the evolution of the 
black hole soft X-ray transients (SXTs) in these articles. On the other hand, 
as acknowledged in them, the observed number of SXTs requires that many, 
if not most, stars of ZAMS masses in the interval of ~ 20 — 35 M© evolve into 
high-mass black holes, of mass ^GM©. 

Brown, Lee, & Bethe (1999) offered as explanation of the large number of 
SXTs that stars in the mass range ~ 20 — 35 Mq could evolve into high- 
mass black holes as long as their H envelope was lifted off by the companion 
star only after the He core burning was complete (Case C mass transfer). If 
it were lifted off earlier, too small an Fe core evolved to make a high-mass 
black hole. In fact, Timmes, Woosley, & Weaver (1996) had found that Fe 
cores of Type lb supernovae were systematically of lower mass than those of 
Type II, basically the same concept. The Brown, Lee, & Bethe scenario had 
been outlined by Brown, Weingartner, & Wijers (1996). Their results were 
not generally believed, however, because the calculation of Woosley, Langer, 
& Weaver (1995) on which they relied employed mass loss rates that were too 
large by a factor of 2 to 3 for the naked He stars. Thus, the present work using 
lower, more correct, rates is necessary. 

Wellstein & Langer (2000) recalculated the evolution of naked He cores to 
their carbon-oxygen cores with wind loss rates lowered by a factor of 2, and 
in some cases by larger factors, finding that because of feedbacks, the total 
wind loss scaled more slowly than linearly with the multiplayer on the wind 
loss rate. In work that we describe here, Alexander Heger, using the Woosley 
& Weaver computer program Kepler, evolved the CO cores further to their Fe 
cores. These arc not sufficiently massive to end up as high-mass black holes; 
in fact, the situation is nearly the same as found by Brown, Weingartner & 
Wijers (1996) where too large He star mass loss rates had been used. Our 
main objective in this article is to explain why the masses of the final Fe 
cores depend only insensitively on the He core mass loss rates. We interpret 
the Fe core mass as a good indicator of the final compact core mass because 
binding energy corrections and mass fallback very nearly compensate each 
other (Brown, Weingartner, & Wijers 1996). 
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The plan of this paper is as follows: 

In Sec. 2 we discuss the improved measurements of wind losses as deduced 
from Wolf-Rayet stars. 

In Sec. 3 we show that the large differences between "naked" and "clothed" 
He star burning can be understood in terms of the different ways in which 
carbon burns in the two cases. 

In Sec. 4 we deal with a possible evolution of the ~ 10 Mq black hole in Cyg 
X-1, suggesting that there may be an analogy with the evolution of transient 
sources; namely that the He core must be kept covered by H during most of 
its burning. This requires a very massive progenitor, a WNL star. 

In Sec. 5 we give simple concluding remarks. 



2 Wind Loss Rates in WNE's (naked He stars) 



It is the mass loss rate proposed by Langer (1989a) for hydrogen-free Wolf- 
Rayet stars, which replaced the previously used assumption of a constant Wolf- 
Rayet mass loss rate in many evolutionary calculations (Schaller et al. 1992; 
cf. WLW 1995), which leads to the small final masses of stars above ~ 35 Mq. 
This semi-empirical rate has been criticised as too high, one argument being 
the existence of massive black hole binaries. Indeed, recently studies of Wolf- 
Rayet stars show that originally measured wind losses have to be corrected 
downwards by a factor of 2... 3 to account for their "dumpiness" (Hamann 
& Koesterke 1998). This is supported by polarisation measurements of the 
Thomson scattering, which depend linearly on the wind density (St. -Louis 
et al. 1993, Moffat & Robert 1994), and is in approximate agreement with 
the rates that would be deduced from the observed rate of increase in orbital 
periods for spherical mass loss 



- = ^ (1) 

in Wolf-Rayet binaries. In V444 Cygni P = 0.202 ± 0.018 s yr~^ was obtained 
by Khaliullin et al. (1984) and Mwr = 9.3 ± 0.5 M© by Marchenko et al. 
(1994), resulting in 



MdynC^lxlO-^Meyr-K (2) 
This is to be compared with the 
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M = 0.75 X 10"^ Mq yr 



-1 



(3) 



obtained by St. -Louis et al. (1993) from the polarisation measurements. In 
later work Moffat & Robert (1994) arrive at a mean of (0.7±0.1) x 10"^ Mq yr'^ 
For WNE (defined here as a Wolf-Rayet phase with vanishing surface hydrogen 
abundance) stars, WLW (1993) used a He wind rate of 



which would give a mass loss rate of 1.6 x Mq yr~^ for a 9.3 WR, a 
factor of 1.6 to 2.1 greater than eqs. (2) and (3). Given the errors in measure- 
ment, we feel safe in saying that the true wind loss rate is not less than 1/3 
to 1/2 the WLW value. Furthermore, we should mention that absolute value 
slope of the mass loss relation Eq. (4) agrees well with the latest empirical 
mass loss rate determinations of Nugis & Lamers (2000). 

However, decreasing the rate to 1/3 does not change our conclusions, as shown 

by recent calculations with such a reduced mass loss rates by Wellstein & 
Langcr (1999), by Fryer ct al. (2001), and by the calculations of Heger de- 
scribed later in this paper. These show that stars with initial masses below 
60 M0 are unlikely to produce black holes more massive than 10 M©, from 
naked hehum stars, although high mass black holes (~ 5 — 10 Mq) may result 
from the ZAMS mass range ~ 20 — 35 Mq employed by Brown, Lee, & Bethe 
(1999) provided Case C mass transfer occurs in the binary. The relative in- 
sensitivity of the final Fe core masses to the wind-loss rate depends somewhat 
intricately on the way carbon burns in a naked He core, as we describe in the 
next section. 

Before we go on to the next section we discuss, however, the possible evolution 
of Cyg X-1, a continuously shining black hole binary where the black hole mass 
is observed to be ^ 10 Mq. Even though massive stars are so prone to lose their 
hydrogen envelopes and then suff'cr extreme wind mass loss as hydrogen-free 
Wolf-Rayet stars, Langer (1987) pointed out that the most massive stars may 
in fact avoid this kind of evolution. The reason for doing so is that these stars, 
upon core hydrogen exhaustion, produce an extended region of intermediate 
hydrogen abundance in between the helium core and the hydrogen-rich enve- 
lope. While the latter is supposed to be lost quickly in a so called Luminous 
Blue Variable (LBV) stage (e.g., Stothers 2000), this intermediate region may 
be sufficiently massive that the star does not manage to blow it away during 
the ensuing so called WNL stage (defined here as a Wolf-Rayet phase with 
non- vanishing surface hydrogen abundance). 

An example for the resulting evolutionary sequence O star LBV WNL 
^ SN is the 100 Mq sequence of Langer & El Eid (1986). Also WLW (1993) 
showed that the final masses of stars above ~ 80 Mq may increase again for 



M = 5 X 1^-^{Mwr/ MqY'^ Mq yr 



-1 



(4) 
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larger initial masses, while stars in the range ~ 35...6OM0 obey the opposite 
trend. Even using the high mass loss rate, their 85 Mq WRB star evolved a 
final He core of 9.71 Mq. We will estimate in Section 4 that this can be brought 
up to ~ 16 Mq with the lower, more correct, mass loss rates. The reason that 
the He core of the WLW 85 Mq star was so much more massive than that of 
their 60 Mq star, was not only the higher mass of the progenitor but the fact 
that the WNL stage of the 85 Mq comprised ~ 40% of the Wolf-Rayet stage, 
whereas that of the 60 Mq star was ~ 25%. During the WNL stage the star is 
covered by hydrogen, so that the mass loss rate is smaller than it would have 
been for a naked He star, and moreover the helium core mass increases with 
time. 

Although the evolution of star as massive as ~ 85 Mq is very uncertain, we 
believe that the only possibility for them to end up in 10 Mq black holes is 
an extended WNL stage, lasting a major fraction of the core helium burning. 
As we discuss below, these stars would have an active hydrogen burning shell 
source, whereas the existence of such a stage in a single massive star would 
probably not affect its evolution substantially, an expansion of the star due to 
hydrogen shell burning does give the possibility of common envelope evolution 
with a companion. In other words, the most massive stars are covered by H 
during much of their core He burning time, and this may enable them to end 
up as high-mass black holes. 



3 "Naked" vs "Clothed" Helium Core Burning 

3.1 Central Carbon Abundances vs Convective Carbon Burning 

As noted earlier, Timmes et al. (1996) and Brown, Weingartner, & Wijers 
(1996) showed that compact cores which followed from "naked" helium stars 
were substantially less massive than followed from those "clothed" during he- 
lium burning with hydrogen envelopes. More recently, this has been explored 
in detail by Wellstein & Langer (1999) and Fryer et al. (2001). An important 
property of the stellar core that is significantly influenced by this difference is 
the carbon abundance in the core at central helium depletion. If its abundance 
is high enough, central carbon burning become exothermic enough to signifi- 
cantly overcome the neutrino losses and burn in a convective core. Otherwise 
the first convective burning of carbon can only appear in a shell source. The 
carbon abundance also has significant effect on the location, duration, and 
extent of the carbon burning shells. For example, as the carbon abundance 
goes to zero, no convective shell burning occurs either. The lower the car- 
bon abundance, the further out the first shell forms. The location of the last 
carbon-burning shell sets the size of the carbon-free core which determines 
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the final evolution of the star and the sizes of the iron core (possible direct 
black hole formation) and the silicon core (possible black hole formation by 
fall-back) . 

There are two major factors that contribute to the resulting carbon abundance 
after core helium depletion. The first is the dependence on the mass of the 
helium core and results chiefly from the different behavior of the ^^C(q;,7)^^0 
reaction, as described clearly by Weaver & Woosley (1993). Since this is im- 
portant to our development, we repeat the argument, supplemented by more 
recent developments. 

Carbon is formed in the triple a-process; since this is a three-body process, 
it depends on the density as p^. Carbon is removed, when possible, by the 
^^C(q!,7)^^0 process, which goes linearly with p. Now, for massive helium 
cores the central density scales roughly as M"^/^. With increasing helium core 
mass this implies that carbon formation, which goes as p^, will be cut down, 
compared with ^^C(a, 7)^^0, which goes as p. Here, we also note that both the 
^^C(q;,7)^^0 reaction and the triple a-process have very similar temperature 
dependence in the temperature regime of central helium burning, 2 x 10* K, 
i.e., the carbon production does not depend on temperature in this regime. 

This means that there is a mass above which the post-helium burning central 
carbon abundance is low enough to skip central carbon burning. Woosley & 
Weaver (1995) found this transition at a mass of 19 Mq. However, this limit 
strongly depends on the still uncertain ^^C(q;, 7)^^0 reaction rate and the 
stellar evolution model, in particular the prescription used for convection or 
mixing processes in general, and also mass loss. On this last issue we elaborate 
later in more detail. 

A bare helium star behaves differently from a clothed one in that its helium 
core does not grow in mass due to hydrogen shell burning, as is the case in 
clothed stars, but rather shrinks due to mass loss from the surface - the whole 
star is only a helium core, a Wolf-Rayet star. Generally, the convective core 
tends to comprise a larger fraction of the helium core, by mass, as the central 
helium abundance decreases, due to the density-dependence outlined above. 
This statement depends, again, somewhat on the description of semiconvec- 
tion, but probably applies to most models with either Schwarzschild convection 
(e.g., Schallcr ct al. 1992), overshooting or fast scmiconvcction (e.g., Woosley & 
Weaver 1995), or rotationally induced mixing (Heger, Langer, Woosley 2000). 
Additionally, the mass fraction of the convective core also increases with the 
mass of helium core. However, bare helium cores (early-type Wolf-Rayet stars) 
experience mass loss rates that are sufficiently strong that the convective core 
actually tends to shrink in the run of its evolution, in particular also towards 
the end of helium burning, rather than grow, as it does in a clothed star. This 
is despite the fact that it still comprises an increasingly larger fraction of the 
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remaining total helium core. The important point is that a growth in mass of 
the convective core injects new helium into this convective core, while when 
the mass is constant or decreases, this injection does not occur. 

As the triple alpha reaction depends on the third power of the helium mass 
fraction it loses against the ^^C(q;,7)^^0 reaction toward the end of central 
helium burning, i.e., carbon is mostly burned rather then produced toward the 
end of central helium burning. That switch typically appears at a central he- 
hum mass fraction of around 10-20%. Most importantly, as can be seen from 
the central carbon abundances at the end of He burning in Table 1, the He 
fraction is too low to burn the final carbon, which remains at a central abun- 
dance ^ 0.30 for all mass stars. This is roughly double the carbon abundance 
necessary for convective carbon burning. It is seen in Table 2 that the central 
carbon abundance in naked He stars goes down only slowly with decreasing 
mass loss rate, although by the time the mass loss is reduced by a factor of 
6 the convective burning stops. In the clothed stars the growth of the core 
and its accompanied injection of helium after this time thus leads to a further 
decrease of carbon as compared to the bare helium cores that do not have this 
additional supply of helium. 

In Fig. 1 we show calculations carried out by Tom Weaver (priv. com. 1995) 
using the KEPLER code (Weaver, Zimmerman, & Woosley 1978; Woosley & 
Weaver 1995) on the central carbon abundance after core helium depletion for 
"clothed" (single) stars without (any) mass loss as function of ZAMS mass. 
The rapid drop in Cc at ZAMS mass Mzams ~ 20 Mq causes the disappear- 
ance of convective carbon burning. 



3.2 Iron Core Masses 

It is generally known that the central entropy rises with the mass of the star; 
so that the more entropy is carried off, the less massive will be the final Fe core. 
This statement can be made more quantitatively by using the Bethe, Brown, 
Applegate, & Lattimer (1979) conclusion that the entropy per nucleon in the 
Fe core is ^ 1, in units of the Boltzmann constant. For a ZAMS 20 Mq, the 
central entropy in our Fe core is 0.76 per nucleon when the He core is burned 
clothed, 0.775 when burned naked; i.e., essentially the same. Thus, the total 
entropy of the Fe core is ~ considering some gradient towards higher 
entropy at the outer layers of the core, where A^jv is the number of nucleons in 
it. In this way one can see directly that the entropy that is carried off during 
the evolution of the Fe core will diminish its mass. 

In Fig. 2 we show iron core masses at the time of iron core collapse for a finely 
spaced grid of stellar masses (Heger, Woosley, Martinez- Pinedo, & Langanke 
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2001). Filled circles and crosses correspond to the core masses of "clothed" 
single stars. The circles were calculated by Tom Weaver (1995, priv. com.) 
using the same physics as in Woosley & Weaver (1995), whereas the crosses 
employ the improved weak rates by Langanke & Martinez-Pmedo (2000) for 
electron capture and beta decay. Although the latter are much smaller than 
the rates by Fuller, Fowler, &: Newman (1985) used by Woosley & Weaver 

(1995) , the final Fe core masses are not much changed, for reasons described 
by Heger, Woosley, Martinez-Pinedo, & Langanke (2001). 

For "clothed" single stars, one sees an anticorrelation between the peak in Fe 
core masses at ZAMS mass ~ 23 Mq and the minimum in Cc at ZAMS mass 
~ 21 Mq. In other words, the peak around 23 Mq in ZAMS masses occurs just 
where the central carbon abundance at the end of helium core burning is at 
its minimum of ~ 10%, too low for convective core burning or the formation 
of a carbon-burning convective shell close to the center of the star. 

When sufficient carbon remains, it will be burned convectively in reactions 
such as ^2C+^2C^24|^g^ 20]^g_^Q, ^^^^ ^ temperature of 70 - 80 keV, sev- 
eral times higher than the ~ 20 keV needed for ^^C(q;, 7)^^0. If the carbon 
abundance is lower, the carbon burning phases are reduced and typically big- 
ger carbon-free cores result for otherwise similar stars. Therefore less energy is 
carried away by neutrinos and the entropy in the core stays higher. More mas- 
sive silicon and iron cores can form in this case (Boyes, Heger, and Woosley 
2001). This is clearly seen in Figures 1 & 2. 

Although the central carbon abundance goes up again for ZAMS mass ~ 
25 Mq, and the Fe core mass goes down, by such high masses the gravitational 
energy of the stellar envelope is large and it is difficult for the shock energy 
after collapse to blow it off, so we beheve that these will go into high-mass 
black holes, as well as the stars of ZAMS mass 20 — 23 Mq. 

3.3 Formation of High-Mass Black Holes 

Although there is a lot of uncertainty in the literature about the ZAMS masses 
which end up in high mass black holes, Bethe & Brown (1999) argued that 
this should occur at the proto compact core mass of 

Mpc~ 1-8 Mq. (5) 

Further uncertainty comes in relating the Fe core mass to the compact core 
mass. It can, however, be seen from Table 3 of Brown, Weingartner, & Wijers 

(1996) that the estimated fallback mass following separation in the supernova 
explosion roughly compensates for the binding energy increase in going from 
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the Fc core to the compact object. Thus, we will assume the compact object 
masses to be the same as those of the Fe cores. In addition to the uncertainty 
in this relation, there is the further uncertainty that the fallback mass, which 
is dependent on the density structure in the layers above the Fe core (Fryer 
1999; Janka 2001), is unknown, since consistent supernova explosions have not 
been achieved. 

None the less, we believe it useful to move ahead with our estimated Mpc 
which indicates that single stars with ZAMS masses ^ 20 Mq go into high- 
mass black holes. Brown, Lee, & Tauris (2001) recently showed that given 

presently accepted wind losses (Schaller et al. 1992) the high mass black holes 
in the transient sources with main sequence companion can only be evolved 
in the region of ZAMS mass ~ 20 Mq, so our choice of Mpc is supported to 
this extent by evolutionary arguments. 

We note briefly that SN 1987A, if it formed a black hole, went into a low mass 
black hole of mass 1.5 Mq < Mbh < 1.8 Mq according to the Brown & Bethe 
(1994) estimates. The black hole would have a low mass in this interval because 
the He envelope was blown off in the delayed explosion. Note from Fig. 2 of 
main sequence masses that there is only a narrow interval from ~ 18 Mq to 
20 Mq in which this could happen given the above interval. 

Schaller et al. (1992) use a ^2c(a, 7)1^0 S-factor of S(300 keV) ~ 100 keV 
barns as compared with the Woosley & Weaver 170 keV barns. They bring 
their central abundance down to 0.16 at the end of core He burning only for 
a ZAMS 25 Mq star, so presumably this would be the mass at which the Fe 
cores would begin to rise rapidly in mass with further evolution. We believe the 
Weaver & Woosley value for the ^^C(q;, 7)^^0 rate to be more correct (Boyes, 
Heger, & Woosley 2001), placing the narrow interval of ZAMS masses from 
which SN 1987A can be evolved correctly. In fact, most recently experiments 
including both El and E2 components have been carried out (Kunz et al. 
2001), obtaining Sf^° = (165 ± 50) keV barns. 



4 Possible Evolution of Cyg X-1 

4-1 Space Velocity vs the Mass Loss in the Formation of Cyg X-1 

According to the calculations of Wellstein & Langer (1999) and Fryer et al. 

(2001), it is unlikely that a black hole as massive as the 10 Mq core (Herrero 
et al. 1995) in Cyg X-1 can be evolved from a naked He star (see Table 1 and 
Fig. 2). Although the calculations of naked He stars were extended up to only 
ZAMS 60 Mq stars, the He winds scale with the 2.6 power of the mass, so 



10 



that higher mass He stars would be expected to behave in the same way. In 
fact, Wellstein & Langer (1999) and Fryer et al. (2001) were unable to evolve 
a 10 Mq black hole, whereas the high space velocity of Cyg X-1, 49 ± 14 km 
s^^ (Kaper et al. 1999) indicates substantial mass loss in a Blaauw-Boersma 
kick (Blaauw 1961; Boersma 1961) in the black hole formation, increasing 
the necessary He core mass. We now estimate this mass loss, following the 
development of Nelemans et al. (1999). It is reasonable to evolve Cyg X-1 
analogously to the soft X-ray black hole transient sources, the difference being 
in the high mass companion star which presumably makes Cyg X-1 shine 
continuously. 

The high space velocity of Cyg X-1 can be explained by mass loss in the black 
hole formation, because this loss from the black hole is somewhat off from 
the center of gravity of the system. Prom Nelemans et al. (1999), the runaway 
velocities from symmetric SNe (See Appendix for the derivation) 

Vsys = 213 X AM X m X P'^Il^c x {^bh + m)-^'^ km s"^ (6) 

where masses are in M0, P in days, and AM = Mue — Mbh with Mne the 
He star mass of the Black hole progenitor. Here Pre-drc is the re-circularized 
orbital period after Blaauw-Boersma kick, and we assume no orbital evolution 
{Pre-circ = Pobs) since the beginning of the mass transfer phase, and neglect 
small eccentricity as in Nelemans et al. (1999). By putting in average values, 

Vsys = 8.32 X AM km s~^. (7) 

In order to obtain the observed velocity, we need AM ~ 5.9 M©, indicating 
the black hole progenitor mass to be Mhc ~ 16 M©. Note that this AM and 
Mhc are in the same ballpark as those estimated by Nelemans et al. for Nova 
Scorpii 1994(GRO J1655-40). This suggests that the evolution of Cyg X-1 is 
similar to that of Nova Scorpii 1994, except that the progenitor mass of the 
black hole must be much higher in Cyg X-1 and the companion is an 0-star 
rather than an F-star. 

The current orbital separation of Cyg X-1 is 



"now 



Before the explosion, the orbital separation of the black hole progenitor and 
companion star was 



"•now no 7-> /r>\ 

OpreSN ~ :j— ~ 33 Rq (9) 

i + 62 
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where 62 = AM/ (Mbh + Mq) ~ 0.21 is the eccentricity right after the super- 
nova explosion. For the derivation, see Appendix. This is only slightly larger 
than the 30 Rq estimated by Bethe & Brown (1999) for the minimum initial 
separation of the 0-star progenitors. If the progenitors were initially closer 
together they would merge already at this stage of evolution. 

As noted by Bethe & Brown (1999), the two 0-star progenitors must have 
initially had a separation of nearly 33 Rq. Substantial mass loss, especially by 
the WR star would have been expected to widen the orbit substantially 



where Mi and Mf are the initial and 
the orbit to the 33 Rq we find before 
there must have been some period of 



(10) 

final system masses. In order to tighten 
the massive star goes into a black hole, 
nonconservative mass transfer. 



4.2 Extended WNL Stage of Black Hole Progenitor 



Langer (1989b) finds that towards the end of core hydrogen burning the radius 
of a massive star may be much larger than its radius on the zero age main 
sequence. For example, a 100 star increases its radius from 13 Rq at H 
ignition by a factor ~ 4 to 53 it!© at core H exhaustion. In the evolutionary 
phase between central H- and He-burning massive stars develop an interme- 
diate fully convective zone just above the H-burning shell which reaches its 
maximal spatial extent and ~ 18 Mq in mass for the 100 Mq star (Langer 
1987). There may be further expansion in a hydrogen shell burning stage un- 
less large mass loss sets in. But the latter may begin only later in the LBV 
stage. In any case, there will be a period before or in the LBV stage where 
the hydrogen is being propelled outwards, but does not yet have enough ve- 
locity to escape. A companion to the very massive star will, if in this region, 
couple hydrodynamically to the hydrogen in the manner explicitly worked out 
by Bethe & Brown (1999), and furnish energy to it from a drop in its gravi- 
tational binding; i.e., through some mechanism resembling common envelope 
evolution. We cannot be more explicit, because of the uncertainty in the post 
main sequence evolution of the very massive stars, but some nonconservative 
mass transfer seems to be necessary in the evolution of a binary such as Cyg 
X-1. We admit to being in somewhat of a predicament, however, since we do 
not want a common envelope evolution which removes the hydrogen before 
the helium burning is well underway. 

It is clear that the black hole in Cyg X-1 cannot have evolved from a naked He 
star. An important clue to its possible evolution is offered by the evolution of 
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WLW (1993) of a ZAMS 85 M© star, where the WR stage began with a WNL 
stage which took ~ 40% of the total WR time. (Their 60 Mq star also had a 
WNL stage, but it took only ~ 25% of the WR time.) During this time the 
WR is covered by hydrogen and the attendant mass loss rate is much lower 
than that for a naked He star. 

Using the times for each WR stage as WLW 1995, we calculate in Table 3 
the masses for ZAMS 60 Mq and 85 Mq stars. Of course there is feedback on 
these masses from the change in times of each stage with altered mass loss 
rate. In order to take the feedback into account, we use the ratios of Mco cores 
calculated by Fryer et al. (2001) for a ZAMS 60 Mq star for various mass loss 
rates to the original WLW (1995) one as plotted in Fig 3. 

We note from Table 3 that a ZAMS 60 star with WR stage mass loss 
rate reduced by 1/3 loses 1.1 Mq during its WNL stage, at an average rate of 
1.5 X 10^^ A/q yr^^. This is an order of magnitude less than it would lose as a 
WNE (bare) He star, and explain why the final He star mass can be as large 
as 13.5 Mq. The WNL mass loss rate is ~ 2.25 times the main sequence mass 
rate for a ZAMS 60 star of Castor, Abbott, & Klein (1975). 

In Fig. 3 we plot our ratios from Table 3 for He star masses alongside the Fryer 
et al. (2001) ratios of CO core masses. The divergence in our favored range 
of 1/2 — 1/3 is not large, so we beheve our procedure in calculating masses is 
justified. We see from Table 3 that we need a mass loss rate of ~ 0.4 times 
the WLW (1995) one in order to lose ~ QMq in the explosion and to be left 
with an ~ IOMq black hole. 

4.3 Discussion on Cyg X-1 Type Objects 

We want to point out that there is a strong metallicity dependance in our 
model. As the winds of WNL stars are likely radiation driven (Hamann et al. 

2000) , we can expect them to weaken with smaller metallicity. Thus, the lower 
limit of the initial mass range from which on the the final stellar mass can grow 
again will be decreasing. Therefore, we expect more Cyg X-1 type systems at 
lower metallicity. We also would like to point out that there is dependance in 
mixing processes, in particular rotation, which can increase the lower limit of 
the initial mass range. 

Using standard population synthesis, Bethe & Brown (1999) estimated that 
there should be ~ 7 Cyg X-1 hke objects in the Galaxy, where we see only 

one (leaving out LMC X-1 and LMC X-3). Although the 7-rays from such 
an object easily penetrate the Galactic disc, they would not be distinguished 
from those from gamma ray bursts (Maarten Schmidt, private communication 

2001) . 
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As possible WR progenitor of the Cyg X-1 black hole, we suggest "WR22"; 
the most massive Wolf-Rayet ever weighed (Rauw ct al. 1996). The minimum 
Wolf-Rayet mass is 72 Mq and the mass ratio to 0-star is 2.78, giving an 0-star 
mass of > 26 Mq. The spectrum of the WN7 exhibits a considerable amount of 
hydrogen {Xh ~ 40%, Hamann et al. 1991, Crowther et al. 1995), strong WR 
emission lines and absorption lines that belong to the WR component (Niemela 
1973; Moffat & Seggewiss 1978). Rauw ct al. (1996) identify the companion 
as an 0-star. With ~ 80 day period, the binary is, however, much too wide in 
order to narrow sufficiently in common envelope evolution to produce a Cyg 
X-1 type binary, according to the Bethe & Brown (1999) estimates. 



5 Conclusion 



Our conclusion to all of the above is simple: In order to evolve high-mass black 
holes at solar metallicity, the He core of the massive star must be covered by 
hydrogen at least most of the time, while it burns. 

For SXTs, the stars in the mass range ~ 20 — 35 Mq could evolve into high- 
mass black holes with the condition that their H envelope was hfted off by 
the companion star only after the He core burning was complete. For more 
massive stars like Cyg X-1 type objects, an extended WNL stage could drive 
the formation of high-mass black holes. 
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Appendix. Blaauw-Boersma Kick 



In this Appendix, tlie details of tlie Blaauw-Boersma kick (Blaauw 1961, 
Boersma 1961) are summarized. Before the explosion,^ from Fig. 4 (a), Rhs, 
VHe, and ui are given as 



/ m \ 

^^f" ^ Tf \ '^1 

XMHe + mj 

VHe = — {GmRHef^ 

ai 

After the explosion, we assume that AM = Mh^ — Mbh "was lost from the 
BH progenitor without interacting with the binary system. The momentum 
lost from the binary system should be compensated by the space velocity of 
the new binary system. 



Ap = AM • VHe = {Mbh + m) ■ v^ys (A.2) 
which gives the space velocity of the new binary system (cm. motion of B) 



= 77 \ ^^^e- (A. 3) 

^ Mbh + m' 



Right after the explosion, in the new cm. frame of BH and 0-star binary, by 
defining the semi-major axis d2 and the eccentricity 62, we have 



UJ2 ■ 



m \ f Gm \ ^^"^ 



Mbh + fnJ \ do 



/ m \ 

RbH = TT ] Ctl 

\Mbh + rrij 

VbH = VHe + Vsys = — ■ VHe- (A.4) 

\MBH+mJ 

Since the vbh is greater than that of the circular motion with radius Rbh 



^ Subscript 1 (2) indicates the properties before (right after) the explosion, and 
subscript "now" indicates the current properties observed. 
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/ MHe + m \ 

vbh = X? T — 

\Mbh + mj 



one can conclude that Rbh is at the minimum of the orbital separation, i.e. 



RBH = d2{l-e2). (A.6) 

Note that vbh and Rbh are perpendicular right after the explosion as in 
Fig. 4 (c). 

The orbital angular momentum of the BH can be expressed using 02 and 62 



la = [Gmn''d2{l - el)] ^'^ = [Gmii^RBH{l + 62)] (A.7) 



1/2 



or 



h = MbhRbhVbh = l^aiVBH = A* {GmRBHf^ {^^^-T~~] (^-8) 
where // is the reduced mass 



Mbh m 

A* = T7 ■ (A-9) 

MsH + m 

By equating la = h we have 



62 = — . A.IO 

Mbh + m ^ ' 



We assume that there is no angular momentum lost during the circulariza- 
tion process by tidal locking after SN explosion. Prom angular momentum 
conservation 



Inow^iGrn/i'^dnow) ^ = [Gm/x^d2(l - 62) . (A.H) 
Therefore, from 



1/2 



Mbh + my 
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-■d2{l-el)^ RBHil + e2) 



m 



~Mbh + m 

we have the relation of the orbital separations 



01(1 + 62), (A.12) 



O-now = ai(l + 62). 

Using this relation, one have 



(A.13) 



CUi = (1 + 62)^ Unow 

or the period relations 



(A. 14) 



(A.15) 



Now Eq. (A. 3) can be expressed in terms of observed quanties using 



VHe 



2/3 



27r 

1 + 62X^/2 



G{Mbh + m) 



1/3 



\ Mhc + 'ITT' , 



1 



Gw? 



^now 



1/2 



{Mbh + m)2/3 p, 



1/3 



Finally the Blaauw-Boersma kick velocity becomes 



(A.16) 



AM (27rGY/^AM m 

VHe - 



Mbh + m {^Mbh + mf'^ P, 



1/3 



:213 



/AM\ / m \ f P, 



\ -^3 / , . X -5/3 

now \ / ^^^BH + m\ 



\MqJ\ MqJ Vdays^ 



Ma 



km s-^ . (A.17) 
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Model 


Mms 


■^comp 




Mhc 


Mco 


Mpe 


Mass Transfer 


Is 


60 


34 


0.30 




3.37 


1.35 


B 


2s 


60 


34 


0.32 


4.07 


3.07 


1.50 


A+AB 


5s 


40 


30 


0.33 


3.84 


2.87 


1.49 


A+AB 


7s 


30 


24 


0.34 


3.63 


2.71 


1.46 


A+AB 


10s 


25 


24 


0.36 


3.39 


2.42 


1.49 


A+AB 


17s 


20 


18 


0.36 


3.39 


2.18 


1.56 


B 



Cc '■ Central Carbon abundance 
Masses : in unit of Mq 

Table 1 

Results of Wellstein & Langcr (1999) and Fryer et al. (2001) with reduced WR mass 
loss rate, 1/2 of WLW (1993) value, at the end of central He burning. The Mpe is 
the final Fe core mass. 



Model Mass Loss Rate* Cc tc Mco[Mq\ MfcIMq] 



Isl 


1 


0.35 


3800 yrs 


2.35 


1.32 


ls2 


1/2 


0.30 


1600 yrs 


3.37 


1.35 


ls3 


1/3 


0.27 


700 yrs 


4.76 


1.61 


ls4 


1/4 


0.25 


500 yrs 


5.93 


1.75 


ls6 


1/6 


0.22 


No 


8.53 


1.50 



* Ratio of WR mass loss to that of WLW 1993 

Cc '■ Central Carbon abundance 

Tc '■ Convective carbon core burning time 

Table 2 

Results of Wellstein & Langer (1999) and Fryer et al. (2001) for different mass 
loss rates of a ZAMS 60 star at the end of central He burning. Models "ls#" 
correspond to the binary system (ZAMS 60 and 34 M©) with Case B mass 
transfer. 
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60 M© WRB 


WR Stage Mass Loss Rate 


r (10^ yrs) 


1 1/2 0.4 1/3 1/4 1/6 


pre WR stage 


26.5 26.5 26.5 26.5 26.5 26.5 


1.13 (25%) 
1.64 (36%) 
1.77 (39%) 


post WNL 


22.5 24.2 24.5 24.8 25.1 25.6 


post WNE 


11.5 15.4 16.5 17.5 18.9 20.8 


post WC/WO 


6.65 9.60 10.7 11.6 13.2 15.5 


after feedback* 


6.65 9.53 11.9 13.5 16.8 24.1 






85 Mq WRB 


WR Stage Mass Loss Rate 


r (10^ yrs) 


1 1/2 0.4 1/3 1/4 1/6 


pre WR stage 


45.3 45.3 45.3 45.3 45.3 45.3 


1.40 (41%) 
1.25 (36%) 
0.80 (23%) 


post WNL 


34.3 40.0 40.8 41.6 42.5 43.4 


post WNE 


14.7 20.9 23.0 24.7 27.5 31.2 


post WC/WO 


9.7 14.3 16.1 17.6 20.2 24.0 


after feedback* 


9.7 13.9 17.3 19.6 24.5 35.2 





* scaled by the numerical results of Fryer et al. (Fig. 3) 



Table 3 

Mass-loss-rate-dependences of masses (in Mq) at the different stages of WR stars 
(models 60WRB and 85WRB of WLW93). The WR stage mass loss rates are in 
units of standard rate of WLW93. The same WR stage times of WLW93 are used 
for different mass loss rates. 
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Fig. 1. Central carbon abundance at the end of He core burning for "clothed" (sin- 
gle) stars as function of ZAMS mass. The calculation was carried out by Tom Weaver 
(1995, priv. com.)using the KEPLER code (Weaver, Zimmerman, & Woosley 1978; 
Woosley k Weaver 1995). The rapid drop in Cc at ZAMS mass Mzams ~ 20 Mq 
signals the disappearance of convective carbon burning. 
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Fig. 2. Comparison of the iron core masses resulting from the evolution of "clothed" 
and "naked" He cores. Filled circles and crosses correspond to the core masses of 
"clothed" stars at the time of iron core implosion for a finely spaced grid of stellar 
masses (Heger, Woosley, Martinez-Pinedo, &: Langanke 2001). The circular were 
calculated with the Woosley & Weaver 1995 code, whereas the crosses employ the 
vastly improved Langanke, Martinez-Pinedo (2000) rates for electron capture and 
beta decay. Open circles (square) correspond to the naked He stars in case A+AB 
(B) mass transfer of Fryer et al. (2001), with reduced WR mass loss rate, see Table 1. 
If the assembled core mass is greater than Mpc = 1.8 Mq, where Mpc is the 
proto-compact star mass as defined by Brown & Bethe (1994), there is no stability 
and no bounce; the core collapses into a high mass black hole. M^s = 1-5 Mq 
denotes the maximum mass of neutron star (Brown & Bethe 1994). The mass of 
the heaviest known well-measured pulsar, PSR B1913+16, is also indicated with 
dashed horizontal line (Thorsett &: Chakrabarty 1999). 
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Fig. 3. Mass-loss-rate-dependences of final CO core masses Mco (filled square) 
from the 60 star discussed in Fryer et al. (2001) and final He core masses (filled 
triangles and circles for 60 Mq and 85 Mq, respectively) after WR stage from Tab. 1. 
The masses are scaled by those with standard rate of WLW93. 
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Before Explosion Right After Explosion 




(a) (b) (c) 

Fig. 4. Coordinates (a) before tire explosion in the cm. frame, (b) right after the 
explosion in the original cm. frame, and (c) in the new cm. frame of black hole 
and companion rhght after the explosion. 
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